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ABSTRACT 

We have obtained long-slit observations in the optical and near infrared of 12 
circumnuclear HE regions (CNSFR) in the early type spiral galaxies NGC 2903, 
NGC 3351 and NGC 3504 with the aim of deriving their chemical abundances. Only 
for one of the regions, the [SIII] A 6312 A was detected providing, together with 
the nebular [SIII] lines at AA 9069, 9532 A, a value of the electron temperature of 
Te([SIII])= 8400^1250^. A semi-empirical method for the derivation of abundances in 
the high metallicity regime is presented. 

We obtain abundances which are comparable to those found in high metallicity 
disc HII regions from direct measurements of electron temperatures and consistent 
with solar values within the errors. The region with the highest oxygen abundance 
is R3+R4 in NGC 3504, 12+log(0/H) = 8.85, about 1.5 solar if the solar oxygen 
abundance is set at the value derived by Asplund et al. (2005), 12+log(O/II)0 = 
8.66±0.05. Region R7 in NGC 3351 has the lowest oxygen abundance of the sample, 
about 0.6 times solar. In all the observed CNSFR the 0/H abundance is dominated 
by the 0+/H+ contribution, as is also the case for high metallicity disc HII regions. 
For our observed regions, however, also the S+/S^+ ratio is larger than one, contrary 
to what is found in high metallicity disc HII regions for which, in general, the sulphur 
abundances are dominated by S^+/H+. 

The derived N/0 ratios are in average larger than those found in high metallicity 
disc HII regions and they do not seem to follow the trend of N/0 vs 0/H which marks 
the secondary behaviour of nitrogen. On the other hand, the S/0 ratios span a very 
narrow range between 0.6 and 0.8 of the solar value. 

As compared to high metallicity disc HII regions, CNSFR show values of the O23 
and the N2 parameters whose distributions are shifted to lower and higher values 
respectively, hence, even though their derived oxygen and sulphur abundances are 
similar, higher values would in principle be obtained for the CNSFR if pure empirical 
methods were used to estimate abundances. CNSFR also show lower ionisation pa- 
rameters than their disc counterparts, as derived from the [SII]/[SIII]. Their ionisation 
structure also seems to be different with CNSFR showing radiation field properties 
more similar to HII galaxies than to disc high metallicity HII regions. 
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1 INTRODUCTION 

The inner IKpc) parts of some spiral galaxies show high 
star formation rates and this star formation is frequently 
arranged in a ring or pseudo-ring pattern around their nu- 
clei. This fact seems to correlate with the presence of bars 
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1985 



and, in fact, computer models which simulate the behaviour 
of gas in galactic potentials have shown that nuclear rings 
may appear as a consequence of matter infall owing to res- 
onances present at the bar edges (Combes & Gerin 
Athanassoula||1992 ). 

In general, Circumnuclear Star Forming Regions 
(CNSFR), also referred to as "hotspots", and luminous and 
large disc HII regions are very much alike, but look more 
compact and show higher peak surface brightness ( KennT] 
cutt et al. 1989 ). In many cases they contribute substantially 



to the emission of the entire nuclear region. 

Their large Ha luminosities, typically higher than 10^^ 
erg s~^, point to relatively massive star clusters as their 
ionisation source, which minimizes the uncertainties due to 
small number statistics when applying population synthe- 
sis techniques (see e.g. Cervi ho et al.n2002| ). These regions 
then constitute excellent places to study how star formation 
proceeds in high metallicity, high density circumnuclear en- 
vironments. 

In many cases, CNSFR show emission line spectra simi- 
lar to those of disc HII regions. However, they show a higher 
continuum from background stellar populations as expected 
from their circumnuclear location, often inside 500 pc from 
the galaxy centre. In early type spirals CNSFR are also ex- 
pected to be amongst the highest metallicity regions as cor- 
responds to their position near the galactic bulge. These 
facts taken together make the analysis of these regions com- 
plicated since, in general, their low excitation makes any 
temperature sensitive line too weak to be measured, partic- 
ularly against a strong underlying stellar continuum. In fact, 
in most cases, the [OIII]A 5007 A line, which is typically one 
hundred times more intense than the auroral [OIII]A 4363 
A one, can be barely seen. 

Yet, despite its difficulty, the importance of an accurate 
determination of the abundances of high metallicity HII re- 
gions cannot be overestimated since they constitute most 
of the HII regions in early spiral galaxies (Sa to Sbc) and 



the inner regions of most late type ones (Sc to Sd) (Diaz 
19891 iVila-Costas & Edmundsll992') without which our de- 



scription of the metallicity distribution in galaxies cannot 
be complete. In particular, the effects of the choice of dif- 
ferent calibrations on the derivation of abundance gradients 
can be very important since any abundance profile fit will 
be strongly biased towards data points at the ends of the 
distribution. It should be kept in mind that abundance gra- 
dients are widely used to constrain chemical evolution mod- 
els, histories of star formation over galactic discs or galaxy 
formation scenarios. The question of how high is the highest 
oxygen abundance in the gaseous phase of galaxies is still 
standing and extrapolation of known radial abundance gra- 
dients would point to CNSFR as the most probable sites for 
these high metallicities. 

Accurate measures of elemental abundances of high 
metallicity regions are crucial to obtain reliable calibrations 
of empirical abundance estimators, widely used but poorly 
constrained, whose choice can severely bias results obtained 
for quantities of the highest relevance for the study of galac- 
tic evolution like the luminosity-metallicity (L-Z) relation for 
galaxies. CNSFR are also ideal cases to study the behaviour 
of abundance estimators in the high metallicity regime. 

In section 2 of this paper we describe the sample of re- 
gions observed and the circumnuclear environment of their 



host galaxies. Our observations and the data reduction pro- 
cedures are described in section 3. The results are presented 
in section 4. The method for the derivation of chemical abun- 
dances is described in section 5. Section 6 is devoted to the 
discussion of our results. Finally, section 7 summarizes the 
main conclusions of this work. 



2 SAMPLE SELECTION 

We have obtained moderate resolution observations of 12 
CNSFR in three "hot spot" galaxies: NGC 2903, NGC 3351 
and NGC 3504 whose main properties are given in Table [l] 
The three of them are early barred spirals and show a high 
star formation rate in their nuclear regions. They are quoted 
in the literature as among the spirals with the highest overall 
oxygen abundances ( Zar itsky et al.|1994) |Perez-01ea,1996| ) . 

NGC 2903 is a well studied galaxy. The Paschen a image 
obtained with the Hubble Space Telescope (HST) reveals the 
presence of a nuclear ring- like morphology with an apparent 
diameter of approximately 15 



625 pc ( Alonso-Herrero 
et al.|[200T ). This structure is also seen, though less promi- 



nent, in the Ha observations from Planesas et al. (1997). A 



large number of stellar clusters are identified on high reso- 
lution infrared images in the K' and H bands, which do not 
coincide spatially with the bright HII regions. A possible in- 
terpretation of this is that the stellar clusters are the result 
of the evolution of giant HII regions (e.g. |Alonso-Herrero 



et al.|[200T ). The global star formation rates in the nuclear 



ring, as deri ved from its Ha luminosity is found to be 0.1 



M, 



^ by[Planesas et al.| ( |l997| and 0.7 



M 



from 



(Alonso-Herrero et al. 2001| . From CO emission observa- 



tions, Planesas et al. (1997) derive a mass of molecular gas 
(H2) of 1.8 X 10 M0 inside a circle 1 Kpc in diameter. 

NGC 3351 is another well known "hot spot" galaxy 
(Sersic & Pastoriza 1967). Early detailed studies of its 



nuclear regions (Alloin & Nieto 1982) concluded that 



NGC 3351 harbours high-mass circumnuclear star forma- 
tion. In fact, the star formation rate per unit area in the 
nuclear region is significantly increased over that observed 



in the disc (Devereux et al.||1992). Elmegreen et al. (1997) 



from near infrared photometry in the J and K bands derive 
a circumnuclear star formation rate of 0.38 yr~^ 



Plane- 



sas et al. (1997), from the Ha emission, derive a total star 
formation rate for the circumnuclear region of 0.24 yr~^ 
and a mass of molecular gas of 3.5 x 10^ Mq inside a circle 
of 1.4 Kpc in diameter, from CO emission observations. A 
recent kinematical study of the CNSFR has been presented 
in iHagele et al.| ( |2007| . 

NGC 3504 is the brightest galaxy in the optically se- 
lected catalogue of starburst galactic nuclei in l Balzano| 
(1983). It forms a pair with NGC 3512. Various studies at 
different wavelengths confirmed that it harbours a very in- 
tense nuclear starburst ( Devereux|1989 Puxley et al.|1990 ). 
Infrared observations in the J and K bands reveal a ring 
with five discrete clumps of st ar formation with colou rs in- 
dicating ages of about 10^ yr ( Elmegreen et al.||l997 ). The 
Ha emission from Plan esas et al.| (1997) traces a compact 
ring structure with a radius of 2 '(200 pc) around the nu- 
cleus where four separated HII regions can be identified. 
From the Ha emission, these authors derive a global star for- 
mation rate for the circumnuclear region of 0.62 M0yr~^, 
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Table 1. The galaxy sample 



Property 



NGC 2903 NGC 3351 NGC 3504 



R. A. (2000)^ 
Dec (2000)" 
Morph. Type 
Distance (Mpc)^ 
pc/ 
Bt (mag)« 
E(B-V),,^(mag) " 



09 32 10.1 
+21 30 03 
SBbc 
8.6 
42 
9.7 
0.031 



10 43 57.7 
+11 42 14 
SBb 
10 
50 
10.5 
0.028 



11 03 11.2 
+27 58 21 
SABab 
20 
100 
11.8 
0.027 



" Ide Vaucouleurs et al. (|l99l|) 


^ NUC 2903: 1 


jBottinelli et al.7l984') 


NGC 3351: 


^aham & et al. ( 199J7| 


NGC 3504:* 


Kenney et al.|(jl992|) 



while from their CO observations, a molecular gas mass of 
21 X 10^ inside a circle 2.7 Kpc in diameter is derived. 



3 OBSERVATIONS AND DATA REDUCTION 

Our spectrophotometric observations were obtained with 
the 4.2m William Herschel Telescope at the Roque de los 
Muchachos Observatory, in 2001 January 26, using the ISIS 
double spectrograph, with the EEV12 and TEK4 detectors 
in the blue and red arm respectively. The incoming light 
was split by the dichroic at A7500 A . Gratings R300B in 
the blue arm and R600R in the red arm were used, cover- 
ing 3400 A in the blue (A3650 to A7000) and 800 A in the 
near IR (A8850 to A9650) and yielding spectral dispersions 
of 1.73 A pixel" ^ in the blue arm and 0.79 A pixel" ^ in 
the red arm. With a slit width of lV05, spectral resolutions 
of ^-^2.0 A and 1.5 A FWHM in the blue and red arms re- 
spectively were attained. This is an optimal configuration 
which allows the simultaneous observation of a given region 
in both frames in a single exposure. 

The nominal spatial sampling is 0^.^4 pixel" ^ in each 
frame and the average seeing for this night was '-^1^.^2. A 
journal of the observations is given in Table [2] 

Two, three and one slit positions in NGC 2903, NGC 
3551 and NGC 3504 respectively were chosen to observe a to- 
tal of 12 CNSFR. Figure [l] shows the different slit positions 
superimposed on images obtained with the HST WFPC2 
camera and taken from the HST archive. Some character- 
sties of the observed regions, as given by Planesas et al. 
(1997), from where the identification numbers have also been 
taken, are listed in Table [3] 

The data were reduced using the IRAl|^ package fol- 
lowing standard methods. The two-dimensional wavelength 
calibration was accurate to 1 A in all cases by means of Cu, 
Ne and Ar calibration lamps. The two-dimensional frames 
were flux calibrated using four spectroscopic standard stars: 
Feige 34, BD26+2606, HZ44 and HD84937, observed before 
and after each programme object with a 3^^ width slit. For 
two of the standard stars: Feige 34 and HZ44, the fluxes have 



^ IRAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatories, which is oper- 
ated by the Association of Universities for Research in Astronomy, 
Inc. (AURA) under cooperative agreement with the National Sci- 
ence Foundation (NSF). 



been obtained fro m the mos t updated version of the origi- 
nal Oke's spectra ( |Oke||l99Q| ) and cover the 3200 to 9200 A 
range. Data between 9200 and 9650 A have been obtained 
from stellar atmosphere models. For the other two stars: 
BD26-h2606 and HD 84937, the fluxes have been taken from 
Oke Gunn] ( |1983| ) that cover the whole spectral range. 
The agreement between the individual calibration curves 
was better than 5% in all cases and a weighted mean cali- 
bration curve was derived. The spectra were previously cor- 
rected for atmospheric extinction using a mean extinction 
curve applicable to La Palma observing site. 

Regarding background subtraction, the high spectral 
dispersion used in the near infrared allowed the almost com- 
plete elimination of the night-sky OH emission lines and, in 
fact, the observed A9532/A9069 ratio is close to the theo- 
retical value of 2.44 in most cases. Telluric absorptions have 
been removed from the spectra of the regions by dividing 
by the relatively featureless continuum of a subdwarf star 
observed in the same night. 



4 RESULTS 

Figure [2] shows the spatial distribution of the Ha flux along 
the slit for the six different positions observed in the sample, 
a single one in the case of NGC 3504, two in NGC 2903 and 
three in NGC 3351. The regions that were extracted into 1-D 
spectra are delimited by arrows. The spectra corresponding 
to each of the identified regions are shown in Figures |3] |4] 
and|5]for NGC 2903, NGC 3351 and NGC 3504 respectively. 



4.1 Underlying population 

The presence of underlying Balmer stellar absorptions is 
clearly evident in the blue spectra of the observed regions 
(see Figures [3] [4] and [5]) and complicates the measurements. 
A two-component - emission and absorption - gaussian fit 
was performed in order to correct the Balmer emission lines 
for this effect. An example of this procedure can be seen 
in Figure [6] The equivalent widths (in amstrongs) of the 
gaussian absorption components resulting from the fits are 
given in Table 4 together with the ratio between the line 
flux measured after subtraction of the absorption compo- 
nent and the line flux measured without any correction and 
using a pseudo-continuum placed at the bottom of the line. 
This factor provides a value for the final correction to the 
measured fluxes in terms of each line flux. In regions R4 of 
NGC 2903 and R5 of NGC 3351 the hne is seen only 
in absorption. In region R5 of NGC 3351 also H7 is seen 
only in absorption. No fitting was performed for these lines, 
hence no correction is listed for them in Table 4. In the 
case of the Hel and Paschen no prominent absorption line 
wings are observed that allow the fitting of an absorption 
component as it was done in the case of the Balmer lines. 
These lines were measured with respect to a local contin- 
uum placed at the their base, which partially corrects by 
underlying absorption. 



4.2 Line intensity measurements 

Emission line fluxes were measured on the extracted spectra 
using the IRAF SPLOT software package, by integrating 
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Table 2. Journal of Observations 



Galaxy 


Spectral range 


Grating 


Disp. 


Spatial resolution 


PA 


Exposure Time 




(A) 




(Apx-i) 


rpx-i) 


n 


(sec) 


NGC 2903 


3650-7000 


R300B 


1.73 


0.4 


105 


2 X 1800 


NGC 2903 


8850-9650 


R600R 


0.79 


0.4 


105 


2 X 1800 


NGC 2903 


3650-7000 


R300B 


1.73 


0.4 


162 


2 X 1800 


NGC 2903 


8850-9650 


R600R 


0.79 


0.4 


162 


2 X 1800 


NGC 3351 


3650-7000 


R300B 


1.73 


0.4 


10 


2 X 1800 


NGC 3351 


8850-9650 


R600R 


0.79 


0.4 


10 


2 X 1800 


NGC 3351 


3650-7000 


R300B 


1.73 


0.4 


38 


2 X 1800 


NGC 3351 


8850-9650 


R600R 


0.79 


0.4 


38 


2 X 1800 


NGC 3351 


3650-7000 


R300B 


1.73 


0.4 


61 


2 X 1800 


NGC 3351 


8850-9650 


R600R 


0.79 


0.4 


61 


2 X 1800 


NGC 3504 


3650-7000 


R300B 


1.73 


0.4 


110 


2 X 1800 


NGC 3504 


8850-9650 


R600R 


0.79 


0.4 


110 


2 X 1800 




NGC3504 




the line intensity over a local fitted continuum. The errors 
in the observed line fluxes have been calculated from the 
expression ai — Gc N^/2[l + EW/(NA)]^/^ where Gi is the 
error in the line flux, Gc represents the standard deviation 
in a box near the measured emission line and stands for 
the error in the continuum placement, N is the number of 
pixels used in the measurement of the line flux, EW is the 
line equivalent width, and A is the wavelength dispersion in 



angstroms per pixel. The first term represents the error in 
the line flux introduced by the uncertainty in the placement 
of the continuum, while the second one scales the signal-to- 



noise in the continuum to the line ( Gonzalez-Delgado et al. 



1994) 



The Balmer emission lines were corrected for the under- 
lying absorption as explained above. Then the logarithmic 
extinction at H/3, c(H/3) was calculated from the Balmer line 
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Table 3. Characteristics of the CNSFR observed 



Galaxy 


Region 


Offsets from centre 


Diameter" 


F(Ra) 






// // 


// 


xlO~^^ erg s~^ cm~^ 


NGC 2903 


R1+R2 


-1.7,+3.5 


4.0 


22.8 




R3 


-1.6,-0.3 


2.0 


2.6 




R4 


-0.3,-3.3 


2.4 


10.1 




R6 


+2.3,+2.4 


2.4 


12.7 


NGC 3351 


Rl 


+0.4,+6.5 


2.4 


12.3 




R2 


-2.6,+2.6 


2.4 


16.0 




R3 


-1.5,-6.5 


2.4 


20.7 




R4 


+0.5,5.4 


2.2 


10.5 




R5 


+2.9,-3.5 


2.4 


4.9 




R6 


+4.1,-0.4 


2.4 


6.1 




R7 


+4.8,+3.6 


1.8 


9.5 


NGC 3504 


R3+R4 


-0.5,-1.7 


1.6 


21.6 



" Size of the circular aperture used to measure fluxes 



Table 4. Equivalent widths of Balmer absorption lines for the observed CNSFR 



Galaxy 


Region 


m (A) 


Corr. factor 


H7 (A) 


Corr. factor 


Rf3 (A) 


Corr. factor 


NGC 2903 


R1+R2 


2.3 


1.15 


3.1 


1.17 


2.9 


1.09 




R3 


2.8 


1.17 


3.6 


1.21 


3.5 


1.09 




R4 


4.5 




3.8 


1.08 


4.8 


1.24 




R6 


3.0 


1.04 


2.6 


1.06 


3.9 


1.18 


NGC 3351 


Rl 


2.1 


1.16 


1.6 


1.14 


3.0 


1.08 




R2 


2.4 


1.30 


1.2 


1.09 


3.2 


1.06 




R3 


2.5 


1.16 


1.2 


1.07 


1.2 


1.03 




R4 


2.1 


1.04 


1.2 


1.08 


2.5 


1.07 




R5 


3.7 




2.5 




3.7 


1.27 




R6 


2.4 


1.32 


2.0 


1.04 


3.5 


1.13 




R7 


2.4 


1.26 


2.1 


1.17 


3.4 


1.09 


NGC 3504 


R3+R4 


4.5 


1.25 


3.6 


1.18 


4.4 


1.10 



decrements assuming the Balmer line theoretical values for 
case B recombination Brocklehurst ( 1971 ) for a temperature 



of 6000 K, as expected for high metallicity regions, and an 
average extinction law |Miller Mathewsl ( 1972| . An exam- 
ple of this procedure is shown in Figure |7| for region R1+R2 
of NGC 2903, where the logarithm of the quotient between 
observed and theoretical Balmer decrements is represented 
against the logarithmic extinction at the Balmer line wave- 
lengths, f(A). The tight relation found for a baseline from 
the Paschen lines to can be taken as evidence of the re- 
liability of our subtraction procedure for this region. This is 
also the case for most regions, although in some cases the 
H7/H;^ ratio lies above or below the reddening line, imply- 
ing that the subtraction procedure is not that good for this 
line. This is not surprising since the presence of the G-band 
and other metal features in the left wing of the line, makes 
the fitting less reliable than for the other Balmer lines. We 
have used all the available Balmer and Paschen-to-Balmer 
ratios, although due to the uncertainties in the emission line 
intensities of the higher order Balmer lines, which are diffi- 
cult to estimate (in fact, in some regions they are seen only 
in absorption) , these ratios have been given a lower weight in 
the fit. This almost amounts to deriving the values of c(H/3) 
from the Ha/H/3 ratio and checking for consistency against 
the values measured for the Paschen lines, which are much 
less affected by underlying absorption due the smaller con- 
tribution to the continuum from main sequence AF stars. 
The errors in c(H/3) have in fact been derived from the mea- 
sured errors in the Ha and H/3 lines, since in any weighted 
average these two lines have by far the largest weight. 



The blue region of the spectrum near the Balmer dis- 
continuity is dominated by absorption lines which cause a 
depression of the continuum and difficult the measurement 
of the [Oil] hues at A 3727 AN one of them however is at 
the actual wavelength of the [Oil] line. This can be seen 
in Figure |8] where we show the spectrum of region R1+R2 
in NGC 2903 compared to that corresponding to a globu- 
lar cluster of M 31, 337-068, of relatively high metallicity 
(Barmby et al. 2000; Mike Beasley, private communication). 
We have measured the line using a local continuum at its 
base as shown in the figure. The different continuum place- 
ments used for computing the error are shown by horizontal 
dotted lines in the figure. 



Once the reddening constant was found, the measured 
line intensities relative to the H/3 line were corrected for in- 
terstellar reddening according to the assumed reddening law. 
The errors in the reddening corrected line intensities have 
been derived by means of error propagation theory. Mea- 
sured and reddening corrected emission line fluxes, together 
with their corresponding errors, are given in Tables |5] |6] and 
|7lfor the observed CNSFR in NGC 2903, NGC 3351 and 
NGC 3504 respectively. Balmer emission lines are corrected 
for underlying absorption. Also given in the tables are the 
assumed reddening law, the H/3 intensity underlying absorp- 
tion and extinction corrected, the H/3 equivalent width, also 
corrected for absorption, and the reddening constant. 



6 A.L Diaz et al 



NGC 2903 PA = 105 



NGC 2903 PA= 162 



4^-15 



Ha linef contmuum 
Ha continuum 
Ha line 





[!ht] 



Pixel number 



Pixel number 



NGC 3351 PA=10 



NGC 3351 PA=38 





1 ' R7 


1 ' 1 


6e-15 


H cell ne+ continuum | 

Ha continuun:! 

— Ha line 


i 




"b 




R 


I 






1 




^ 3e- 15 








i 

n 


1 ■ 




R5 

U 

V.I — . , , . 



' 1 ' 


R2 1 


Ha line+ continuum 


1 


Ha continuum * 




— Ha line 




Rl 

hi 




y 

J.. 





450 

Pixel number 



300 450 
Pixel number 



NGC 3351 PA=61 



Ha line+continuum 
Ha continuum 
a line 




300 450 
Pixel number 



NGC 3504 PA=110 



Ha line+continuum 
a continuum 
Ha line 




300 



Pixel number 



600 



Figure 2. Ho; profiles for the observed slit positions. Each figure includes the name of the galaxy, the P.A. of the slit, the name of the 
observed regions and their X,Y positions (see also Table 2). 



5 CHEMICAL ABUNDANCES 

Electron densities for each observed region have been de- 
rived from the [Sll] A A 6717, 6731 A li ne ratio, following 
standard methods (e.g. OsterbrocklllQSQ ). They were found 
to be, in all cases, < 600 cm~ , higher than those usually 
derived in disk HII regions, but still below the critical value 
for coUisional de-excitation. 

The low excitation of the regions, as evidenced by the 
weakness of the [OIII] A 5007 A line (see left panels of Fig- 
ures [3] [4] and |5| , precludes the detection and measurement 
of the auroral [OIII] A 4363 A necessary for the derivation of 



the electron temperature. It is therefore impossible to obtain 
a direct determination of the oxygen abundances. Empirical 
calibrations have to be used instead. 

Different calibrators for strong emission lines have been 
proposed in the literature for different kinds of objects in- 
volving different chemical elements, among others, the oxy- 
gen abundance parameter R23 = O23 (Pagel et al. 1979), 



the nitrogen N2 parameter ( |DenicolQ et al. 2002) and the 
sulphur abundance parameter S23 (Diaz & Perez-Montero 
2000|). Also, calibrators involving a combination of emi- 



sion lines of two elements have been proposed, such as 
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Figure 3. Extracted blue (left) and red (right) spectra for the observed regions of NGC 2903. From top to bottom: R1+R2, R3, R4 and 
R6. 
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ure 4. Extracted blue (left) and red (right) spectra for the observed regions of NGC 3351. From top to bottom: Rl, R2, R3 and R4. 
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Figure 4. (cont) Extracted blue (left) and red (right) spectra for the observed regions of NGC 3351. From top to bottom: R5, R6 and 
R7. 



[OIII]/[NII] dAlloin et al^|l979t , [ArIII]A 7135 A/ [QUI] A 
5007 A and [SIII]A 9069 A/[OIII] A 5007 A ( |Stasinska|20"06 ). 
Our CNSFR, show very weak lines of [OIII] which are mea- 
sured with large errors. This is taken as evidence for high 
metallicity in these regions. This, in turn, may imply val- 
ues of N/0 larger than solar, due to the chemical evolution 
of the regions themselves and the increasing production of 
secondary nitrogen. Hence, we have considered the use of 
the N2 parameter unreliable for this kind of objects. On the 
other hand, the [SIII] lines are seen to be strong as compared 
to the [OIII] lines (see right panels of Figures [s] [4] and [s]). Re- 



cently, the combination of both the oxygen abundance, O23, 
and sulphur abundance, S23, parameters has been claimed 
to be a good metallicity indicator for high metallicity HII 
regions ( Diaz Perez-Montero|2000 Perez-Montero & Diaz 
2005). From now onwards we will call this parameter SO23, 



and is defined as: 

^3 _ /([5'//]A6716,31) + /([5'///]A9069,9532) 
O23 ~ 



SO, 



/([0//]A3727, 29) + /([0///]A4959, 5007) 
This parameter is similar to the S3O3 proposed by 



Stasinska (2006) but is, at first order, independent of geo- 



metrical (ionization parameter) effects. The amount of avail- 
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Figure 5. Extracted blue (left) and red (right) spectra for regions R3 + R4 of NGC 3504. 
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Figure 6. Examples of the fitting procedure used to correct the Balmer emission line intensities for underlying absorption. 



able data on sulphur emission lines is increasingly growing, 
especially in the high metallicity regime. This makes possi- 
ble for the first time to calibrate the electron temperature 
of [SIII] in terms of the SO23 parameter. To perform this 
calibration we have compiled all the data so far at hand 
with sulphur emission line data both for the auroral and 
nebular lines at A 6312 A and AA 9069,9532 A respectively. 
The sample comprises data on galactic (Garcia- Roj as 2006) 
and extragalactic ( |Diaz et al. 1 19871 IVilchez et al.|1988| |Pas- 



toriza et 



et al. 



Caste 



2004 



1994^ 



aL]|1993| IGarnett et al.| »1997"i Gonzalez-Delgado' 
41 IGonzilez-Delffado et al. 1995; Diaz et al. 2000, 



lanos et al. |2002| | IKennicutt et al. 2003 ; Bresolin et al, 



2005^ HII regions and HII ga laxies (jSkillman Ken- 



nicutt||1993l [Skillman et al.||1994| |Perez-Montero Diaz 



2003; Hagel e et al.| |2006; Hagele G. F. 2007). The data on 
extragalactic HII regions has been further split into low 
and high metalicity HII regions according to the criterion 
oflDiaz & Perez-Monterol (|2000|), i.e. log023 < 0.47 and - 



0.5< log S23 < 0.28 implies oversolar abundances. For all 



the regions the [SIII] electron temperature has been derived 
from the ratio between the auroral and the nebular sulphur 
lines, using a five-level atom program dShaw k Dufour|1995| ) 
and the collisional strengths from [Tayal Gupta (1999), 
through the task TEMDEN as implemented in the IRAF 
package. The calibration is shown in Figure [9] together with 
the quadratic fit to the high metallicity HII region data: 

te{[SIII]) = 0.596 - 0.283/o^^5'O23 + 0.199(/o^^5'O23)^ 

Figure [To] shows the comparison between the electron 
temperatures and ionic sulphur abundances derived from 
measurements of the auroral sulphur line A 6312 A and our 
derived calibration. The S^/H^ and S^^/H^ ionic ratios 
have been derived from: 

12 + log§^= log {'-^^11^^) + 

5.423 + ^ - 0.28 • logt 
12 + /o,|±= /o,(M±f-^) + 

5.8+ ^ -0.22 -logt 
These expressions have been derived by performing ap- 



Abundances of circumnudear star forming regions 11 




Figure 7. Reddening determination for region R1+R2 in NGC 2903. The tight relation found shows the goodness of the correction to 
the Balmer emission hnes by the underlying absorption continuum. 




Figure 8. Blue spectrum of region R1+R2 in NGC 2903 showing the location of the [Oil] AA 3727,29 A hnes. The spectrum below 
corresponds to a metal rich globular cluster in M 31 which can be compared with the underlying stellar population in the region. The way 
in which the [Oil] line has been measured by placing a local pseudo-continuum at its base is shown. The different continuum placementes 
used for the computation of the errors are shown by horizontal lines. 



propriate fittings to the IONIC task results following the 



functional form given in Pagel et al. (1992). The assump- 
tion has been made that Te[SIII] ~Te([SII]) in the observed 
regions. This assumption seems to be justified in view of the 
results presented by Bresolin et al. ( 2005| ). 

We have used the calibration above to derive te([SIII]) 
for our observed CNSFR. In all cases the values of log023 
are inside the range used in performing the calibration, thus 
requiring no extrapolation of the fit. These temperatures, 
in turn, have been used to derive the S^/H^ and S^^/H^ 
ionic ratios. The derived Te and ionic abundances for sul- 
phur are given in Table [S] together with measured values of 
the electron density. Temperature and abundance errors are 
formal errors calculated from the measured line intensity er- 
rors applying error propagation formulae, without assigning 
any error to the temperature calibration itself. 
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Table 7. Reddening corrected emission line intensities for the 
CNSFR in NGC 3504 



Region 






R3+R4 


A (A) 


Line 


fA 


^ X 


T , 


3727 


[Oil] 


0.271 


yoo zt yo 


1 QQC _|_ 1 /I 1 


4102 


Ha 


0.188 


ZOUztO 


oZDzt / 


4340 




0.142 


4touzizo 


r:fi-\ _|_o 


4686 


Hell 


0.045 






4861 


H,3 


0.000 


iUUUzto 


iUUUzto 


4959 


[OIII] 


-0.024 




fil -\-A 

DlZlZ'4: 


5007 


[OIII] 


-0.035 


1 QO-UQ 


1 QQ-l-7' 

ioozt I 


5199 


[NI] 


-0.078 


oyztio 


OZztio 


5876 


Hel 


-0.209 


87±11 


65±8 


6300 


[OI] 


-0.276 


80±7 


54±5 


6312 


[SIII] 


-0.278 






6548 


[Nil] 


-0.311 


818±17 


529±13 


6563 




-0.313 


4423±29 


2850±19 


6584 


[Nil] 


-0.316 


2466±29 


1583±30 


6678 


Hel 


-0.329 


17±4 


11±3 


6717 


[SII] 


-0.334 


652±13 


408±10 


6731 


[SII] 


-0.336 


599±12 


374±9 


8863 


Pll 


-0.546 






9016 


PIO 


-0.557 


34±7 


15±3 


9069 


[SIII] 


-0.561 


327±13 


149±7 


9230 


P9 


-0.572 


95±8 


43±4 


9532 


[SIII] 


-0.592 


802±25 


350±15 


9547 


P8 


-0.593 






c{m 






0.61±0.02 




I(H/3)« 








20.8±0.17 


EW(H^)(A) 






15.2±0.2 


"in units of 10"^^ 


erg s"-*^ 


cm~^ 





Once the sulphur ionic abundances have been derived, 
we have estimated the corresponding oxygen abundances. In 
order to do that, we have assumed that sulphur and oxygen 
electron temperatures follow the relation given by Garnett 
(1992) and confirmed by more recent data (Hagele et al. 
|2006| , and hence we have derived te([OIII]) according to 
the expression: 



te{[OIII]) = 1.205 • te{[SIII]) - 0.205 

We have also assumed that te ([Oil]) ~ te([SIII]). The values 
of Te ( [OIII] ) and the ionic abundances for oxygen are given 
m Table [9] The same comments regarding errors mentioned 
above apply. Finally, we have derived the N^/0^ ratio as- 
suming that te([OII]) te([NII]) ^ te([SIII]). Thcsc values 
are also listed in Tabled 

For one of the observed regions, R1+R2 in NGC 2903, 
the [SIII] A 6312 A has been detected and measured. Al- 
though the line intensity ratios for this region presented in 
Table [5] have been measured on the combined spectrum of 
slit positions 1 and 2, the placement of the continuum for 
the [SIII] 6312 line was very uncertain. Therefore, we have 
performed the temperature measurement on the spectrum 
extracted from slit position 1, which shows the best defined 
continuum. On this spectrum we have measured the inten- 
sity of the [SIII] A 6312 A line with respect to Ha, and 
those of [SIII] AA 9069, 9532 A with respect to P9 A 9329, 
in order to minimize reddening corrections. The region of 
this spectrum around the [SIII] A 6312 A line is shown in 
Figure^^where the horizontal lines show the different place- 
ments of the continuum used to measure the line and cal- 



culate the corresponding errors. The obtained [SIII] nebular 
to auroral line ratio is: 79.4zb49.1 which gives a Te([SIII]) = 
84001^250^5 slightly higher than predicted by the proposed 
fit. This region is represented as a solid black circle in Figure 

El 



6 DISCUSSION 

6.1 Characteristics of the observed CNSFR 

Planesas et al. (1997) provide Ha fluxes (from images) for 
all our observed regions. These values are larger than those 
measured inside our slit by factors between 1.0 and 2.7, 
depending on the size of the region, something to be ex- 
pected given the long-slit nature of our observations. We 
have calculated the Ha luminosities for our regions from our 
observed values, correcting for extinction according to the 
values found from the spectroscopic analysis. The resulting 
values are listed in Table |1Q[ These values are larger than 
the typical ones found for disc HII regions and overlap with 
those measured in HII galaxies. The region with the largest 
Ha luminosity is R3+R4 in NGC 3504, for which a value of 

2.02 X 10^° erg s~-^ is measured. 

We have derived the number of hydrogen ionizing pho- 
tons from the extinction corrected Ha flux as: 



/osQ(//°) = 0.802 X 10^^ (^ 



where F(Ha) is in erg cm~^s~^ and the distance, D, is in 
Mpc. 

The ionization parameter, u, can be estimated from the 
[SII] /[SIII] ratio ( |Dfaz et al.|[lM| as: 

logu = -1.6Slog{[SII]/[SIII]) - 2.99 

and ranges between -3.12 and -3.98 for our observed CNSFR, 
in the low side of what is found in disc HII regions, even in 
the cases of high metallicity (see jDiaz et al.||1991 ). 

From the calculated values of the number of Lyman a 
photons, Q(H°), ionisation parameter and electron density, 
it is possible to derive the size of the emitting regions as well 
as the filling factor (see |Castellanos et al.|2002"] ). The derived 
sizes are between 1.5 arcsec for region R3 in NGC 3351 and 
5.7 arcsec for region R4 in NGC 2903; these values corre- 
spond to linear dimensions between 74 and 234 pc. The de- 
rived filling factors are low: between 6 x 10~^ and 1 x 10~^, 
lower than commonly found in giant HII regions (~ 0.01). 
Sizes in arcsec, filling factors and the corresponding masses 



of ionised hydrogen, M(HII), are given in Table 10 

We have also derived the mass of ionising stars, M*, 
from the calculated numer of hydrogen ionising photons with 



the use of evolutionary models of ionising clusters (Garcia 
Vargas Dfaz'1994"i "Stasinska Leitherer||1996 ) assummg 
that the regions are ionisation bound and that no photons 
are absorbed by dust. A Salpeter IMF with upper and lower 
mass limits of 100 and 0.8 Mq respectively, has been as- 
sumed. According to these models, a relation exists between 
the degree of evolution of the cluster, as represented by its 
}if3 emission line equivalent width and the number of hydro- 



gen ionising photons per unit solar mass ( |Dfaz 1998 Diaz 



et al. 2000). The ionising cluster masses thus derived are 
given in Table 10 and range between 1.1 x 10^ and 4.7 x 



10 M0. The measured H/3 equivalent widths, however, are 
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Figure 9. Empirical calibration of the [SIII] electron temperature as a function of the abundance parameter SO23 defined in the text. 
The solid line represents a quadratic fit to the high metallicity HII region data. References for the data are given in the text. 




Figure 10. Comparison between the electron temperatures (left) and ionic abundances for sulphur (right) for the high metallicity HII 
region sample, as derived from measurements of the auroral sulphur line A 6312 A (abscissa) and from our derived calibration (ordinate). 
The solid line represents the one-to-one relation. 



very low and could be reflecting the contribution by underly- 
ing non-ionising populations. An alternative way to take into 
account the cluster evolution in the derivation of the mass is 
to make use of the existing relation between the ionisation 
parameter and the H/3 equivalent width for ionised regions 
( Hoyos D la^ |2006 |). In that case, the derived masses are 
lower by factors between 1.5 and 15. At any rate, given the 
assumptions of no dust absorption or photon leakage, these 
masses represent lower limits. 



6.2 Metallicity estimates 

The abundances we derive using our Te([SIII]) calibration 
are comparable to those found by Bresolin et al. (2005) for 
their sample of high metallicity HII regions. Most of our 
CNSFR show total oxygen abundances, taken to be 0/H= 
0^/H^ + 0^^/H^, consistent with solar values within the 
errors. The region with the highest oxygen abundance is 
R3+R4 in NGC 3504: 12+log(0/H) = 8.85, about 1.5 solar 
if the solar oxygen abundance is set at the value derived by 
Asplund et al. (2005), 12+log(O/H)0 = 8.66±0.05. Region 



R6 in NGC 3351 has the lowest oxygen abundance of the 
sample, about 0.6 times solar. In all the observed CNSFR the 
0/H abundance is dominated by the 0^/H^ contribution 
with 0.18 < log(0+/0^+) < 0.85. This is also the case for 
high metallicity disc HII regions where these values are even 
higher. For our observed regions, also the S^/S^^ ratios are 
larger than one, which is at odds with the high metallicity 
disc HII egions for which, in general, the sulphur abundances 
are dominated by S^^/H^. 

The fact that both 0/H and S/H seem to be dominated 
by the lower ionisation species, can raise concern about our 
method of abundance derivation, based on the calibration of 



the Te([SIII]). In a recent article [Pilyugin (2007) addresses 
this particular problem. He proposes a calibration of the ra- 
tio of the [Nil] nebular-to-auroral line intensities in terms 
of those of the nebular oxygen lines. Then the [Nil] elec- 
tron temperature, thought to properly characterise the low 
ionisation zone of the nebula, can be obtained. We have ap- 
plied Pilyugin (2007) to our observed CNSFR as well as the 
high metallicity HII sample. Figure ^] shows the derived 
te([NII]) following Pilyugin's method against the te([SIII]) 
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Table 8. Derived electron densities, electron temperatures and ionic abundances from sulphur lines in our observed CNSFR. 



Galaxy 


Region 


Ue (cm~^) 


Te([SIII]) (K) 


12±log 


(S+/H+) 


12±log (S2+/H+) 


12±log[ (S+±S2+)/H+] 


NGC 2903 


R1+R2 


280 ± 90 


5731 ± 188 


6.83 


± 0.04 


6.50 ± 0.05 


7.01 ± 0.05 




R3 


230 ±135 


5737 ± 210 


6.91 


± 0.06 


6.60 ± 0.06 


7.09 ± 0.06 




R4 


270 ±120 


6660 ± 331 


6.68 


± 0.06 


6.21 ± 0.08 


6.81 ± 0.06 




R6 


250 ± 80 


5413 ± 221 


6.97 


± 0.05 


6.64 ± 0.06 


7.14 ± 0.05 


NGC 3351 


Rl 


360 ± 80 


5683 ± 188 


6.74 


± 0.04 


6.75 ± 0.05 


7.05 ± 0.04 




R2 


440 ±110 


5405 ± 210 


6.90 


± 0.05 


6.65 ± 0.06 


7.09 ± 0.05 




R3 


430 ± 70 


5399 ± 210 


6.84 


± 0.04 


6.83 ± 0.05 


7.14 ± 0.05 




R4 


310 ±120 


5390 ± 271 


6.87 


± 0.06 


6.64± 0.07 


7.07 ± 0.07 




R5 


360 ±230 


6177 ± 286 


6.78 


± 0.07 


6.49 ± 0.09 


6.97 ± 0.08 




R6 


360 ±170 


6297 ± 295 


6.49 


± 0.06 


6.20 ± 0.07 


6.67 ± 0.07 




R7 


410 ±100 


5537 ± 220 


6.80 


± 0.05 


6.52 ± 0.07 


6.99 ± 0.06 


NGC 3504 


R3+R4 


370 ± 60 


6288 ± 170 


6.86 


± 0.03 


6.77 ± 0.03 


7.12 ± 0.03 



Table 9. Derived Te and ionic abundances for oxygen and nitrogen in our observed CNSFR. 



Galaxy 


Region 


Te([OIII]) (K) 


12±log (0+/H+) 


12±log (02+/H+) 


12±log[ (0+±02+)/H+] 


log(N+/0+) 


NGC 2903 


R1±R2 


4855 ± 226 


8.55 ± 0.08 


8.11 ± 0.08 


8.69 ± 0.08 


-0.37 ± 0.07 




R3 


4863 ± 253 


8.55 ± 0.09 


8.37 ± 0.09 


8.77 ± 0.09 


-0.36 ± 0.08 




R4 


5975 ± 399 


8.42 ± 0.10 


8.04 ± 0.14 


8.57 ± 0.11 


-0.50 ± 0.08 




R6 


4473 ± 267 


8.59 ± 0.10 


8.33 ± 0.10 


8.79 ± 0.10 


-0.38 ± 0.07 


NGC 3351 


Rl 


4798 ± 227 


8.60 ± 0.08 


8.06 ± 0.08 


8.72 ± 0.08 


-0.49 ± 0.07 




R2 


4463 ± 253 


8.47 ± 0.09 


8.41 ± 0.10 


8.74 ± 0.09 


-0.23 ± 0.07 




R3 


4456 ± 253 


8.64 ± 0.09 


8.17 ± 0.10 


8.77 ± 0.09 


-0.38 ± 0.07 




R4 


4445 ± 327 


8.55 ± 0.11 


8.19 ± 0.13 


8.71 ± 0.12 


-0.37 ± 0.08 




R5 


5393 ± 344 


8.42 ± 0.10 


8.34 ± 0.10 


8.68 ± 0.10 


-0.37 ± 0.10 




R6 


5538 ± 356 


8.34 ± 0.10 


7.62 ± 0.17 


8.42 ± 0.11 


-0.47 ± 0.08 




R7 


4622 ± 266 


8.48 ± 0.09 


8.15 ± 0.10 


8.65 ± 0.10 


-0.40 ± 0.07 


NGC 3504 


R3±R4 


5527 ± 205 


8.79 ± 0.07 


7.94 ± 0.06 


8.85 ± 0.07 


-0.65 ± 0.06 



derived from our calibration. The (red) dashed line shows 
the one-to-one relation while the (black) solid line shows 
the actual fit to all the data. It can be seen that, in what re- 
gards CNSFR, both temperatures are very similar, with the 
[Nil] temperature being, in average, 500 K higher than that 
of [SIII]. This difference is of the same size of the average 
errors in the measured temperatures entering the calibra- 
tions, therefore we can assume that our derived te([SIII]) 
characterises the low ionisation zone at least as well as the 
te([NII]) derived applying Pilyugin's method. 

Concerning relative abundances, with our analysis it is 
possible to derive the relative N/0 value, assumed to be 
equal to the N^/0^ ratio. These values are in all cases larger 
than the solar one (log(N/O)0= -0.88; Asplund et al. 2005) 
by factors between 1.7 (R3+R4 in NGC 3504) and 4.5 (R2 
in NGC 3351) which are amongst the highest observed N/0 
ratios (see e.g. Molla et al.|2006| ). Regarding S, if - given the 
low excitation of the observed regions - the fraction of S^^ 
is assumed to be negligible, the S/0 ratio can be obtained 
as: 

s_ _ + s^^ 

o~ 0+ + 02+ 

The values of log(S/0) span a very narrow range between 
-1.76 and -1.63, that is between 0.6 and 0.8 of the solar value 
(log(S/O)0 = -1.52; [Asplund et al.|[20Q5] ). 



6.3 Comparison with high metallicity HII regions 

The observed CNSFR being of high metallicity show how- 
ever marked differences with respect to high metallicity disc 
HII regions. Even though their derived oxygen and sulphur 



abundances are similar, they show values of the O23 and the 
N2 parameters whose distributions are shifted to lower and 
higher values respectively with respect to the high metallic- 



ity disc sample (Figure 13). Hence, if pure empirical meth- 
ods were used to estimate the oxygen abundances for these 
regions, higher values would in principle be obtained. This 
would seem to be in agreement with the fact that CNSFR, 
when compared to the disc high metallicity regions, show 
the highest [Nil] /[OH] ratios. Figure 14 shows indeed a bi- 
modal distribution of the [Nil] /[OH] ratio in disc and cir- 
cumnudear HII regions. 

A good correlation has been found to exist between 
the [Nil] /[OH] ratio and the N^/0^ ionic abundance ra- 
tio, which in turn can be assumed to trace the N/0 ratio 
( Perez-Montero Diaz|2005 ). This relation is shown in Fig- 
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for the observed CNSFR (black circles) and the high 
metallicity HII region sample (grey circles). The CNSFR 
sample has been enlarged with two circumnudear regions 
of NGC 1097 observed by PhilUps et al.| ([T984| and other 
three observed in NGC 5953 by [Gonzalez Delgado &: Perez] 



( 1996 ), all of them of high metallicity. Also shown are data of 



some CNSFR in two peculiar galaxies: NGC 3310 (Pastoriza 
eraL][l993t and NGC 7714 ( G onzalez-Delgado et al.||1995| ) 



of reported lower metallicity. In all the cases, ionic and total 
abundances have been derived following the same methods 
as in the CNSFR in the present study and described in sec- 
tion 5. In the case of the latter regions, abundances derived 
in this way are larger than derived from direct determina- 
tions of te([OIII]) since our estimated values for this temper- 
ature are systematically lower by about 1200 K. This could 
be due to the fact that our semi-empirical calibration has 
been actually derived for high metallicity regions. However, 



18 A.I. Diaz et al. 




Figure 11. Spectrum of region R1+R2 in NGC 2903 around [SIII] A 6312 A. The way in which the [SIII] hne has been measured is 
shown. The different continuum placements used for the computation of the errors are shown by horizontal lines. 




0.4 0.5 0.6 0.7 O.S 0.9 1 



t/[Sin]) (this work) 

Figure 12. te([NII]) derived from Pilyugin's method against the te([SIII]) derived from the calibration presented in this work. The 
dashed line shows the one-to-one correspondence. 



the values of log(S023) for these regions are between -0.44 
and -0.85 and therefore within the validity range of the cal- 
ibration and we have preferred to use the same method for 
consistency reasons. 

We can see that a very tight correlation exists which al- 
lows to estimate the N/0 ratio from the measured [Nil] and 
[Oil] emission line intensities. In this relation, high metal- 
licity regions and CNSFR seem to follow a sequence of in- 
creasing N/0 ratio. A linear regression fit to the data yields 
the expression: 

log{N/0) = (0.65 ± Om)log{[NII]/[OII]) - (0.79 ± 0.01) 
This relation is shallower than that found in IPerez-Monterol 



& Diaz ( 2005 ) for a sample that did not include high metal- 



licity HII regions. The N/0 ratio for our observed CNSFR 
is shown in Figure [16] against their oxygen abundance to- 
gether with similar data for the high metallicity HII region 
and HII galaxy samples. It can be seen that all the CNSFR 
show similar oxygen abundances, with the mean value being 
lower than that shown by high metallicity disc HII regions, 
but the observed CNSFR show larger N/0 ratios and they 
do not seem to follow the trend of N/0 vs 0/H which marks 
the secondary behaviour of nitrogen. 

The values of the oxygen abundance that we find for 
the CNSFR in NGC 3351 are equal within the errors, with 



a mean value 12+log(0/H) = 8.70 =b 0.10, except for re- 
gion R6 which shows an 0/H abundance lower by a factor 
of about 2. It is worth noting that this is the region with 
the highest 0^/0 ratio: 0.83 and a very low 0^^/H^ value: 
12+log(0^+/H+) = 7.62 it 0.17. The average value for the 
rest of the regions is in agreement with the central abun- 



dance found by Pilyugin et al. (2006) by extrapolating the 
galaxy oxygen abundance gradient, 12+log(0/H) = 8.74 ± 
0.02. This value is somewhat lower than that previously esti- 
mated by Pilyugin et al. ( 2004 ): 8.90 from the same data but 



following a different method of analysis. In that same work, 
the quoted central abundance for NGC 2903 is 12+log(0/H) 
= 8.94. Our results for regions R1+R2, R3 and R6 in this 
galaxy yield an average value of 8.75 ± 0.09 lower than theirs 
by 0.2 dex. Region R4 shows a lower oxygen abundance but 
still consistent within the errors with the average. Values 
of N/0 ratios for the centres of NGC 3351 and NGC 2903 
similar to those found here (-0.33 and -0.35 respectively) are 
quoted by [Pilyugin et al.| ( |2004 ). 



Another difference between the high metallicity circum- 
nuclear and disc regions is related to their average ionisation 
parameter. The left panel of Figure pT] shows the distribution 
of the [SII]/[SIII] ratio for the two samples. The [SII]/[SIII] 
ratio has been shown to be a good ionisation parameter in- 
dicator for moderate to high metallicities ( Diaz et al.|199l] ) 
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Figure 13. Distribution of the empirical abundance parameters O23 (left) and N2 (right) for the observed CNSFR and the sample of 
high metallicity disc HII regions. 
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Figure 14. Distribution of the [Nil] /[Oil] ratio for the observed CNSFR (dark) and the sample of high metallicity disc HII regions 
(light). 








Figure 15. The relation between the [Nil] /[Oil] emission line intensity ratio and the N+/0+ ratio for CNSFR (black cicles) and high 
metalHcity HII regions (grey circles). Included also (downward triangles) are lower metallicity CNSFR in NGC 3310 ( [Pastoriza et al.| 
[19931 and NGC 7714 (Gonzalez-Delgado et al.|1995j . 
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Table 10. General properties of the observed CNSFR. 



Galaxy 


Region 


F(Ha) 


L(Ha) 


Q(hO) 


log u 


Diameter 


e 


M* 


M(HII) 






(erg cm~^ s"-*-) 


(erg s-i) 


(photons s"-*^) 




(arcsec) 




(Mo) 


(Mo) 


NGC 2903 


R1+R2 


3.06E-13 


2.71E+39 


1.98E+51 


-3.69 


4.66 


5.93E-04 


7.39E+05 


1.44E+04 




R3 


3.16E-14 


2.80E+38 


2.05E+50 


-3.65 


1.59 


2.29E-03 


l.llE+05 


1.81E+03 




R4 


2.23E-13 


1.98E+39 


1.45E+51 


-3.98 


5.70 


2.54E-04 


2.11E+06 


1.09E+04 




R6 


1.48E-13 


1.31E+39 


9.60E+50 


-3.66 


3.33 


9.89E-04 


4.55E+05 


7.79E+03 


NGC 3351 


Rl 


2.19E-13 


2.65E+39 


1.94E+51 


-3.12 


1.81 


3.76E-03 


6.40E+05 


1.09E+04 




R2 


1.61E-13 


1.95E+39 


1.42E+51 


-3.55 


2.29 


9.08E-04 


6.80E+05 


6.56E+03 




R3 


2.35E-13 


2.84E+39 


2.08E+51 


-3.12 


1.72 


3.31E-03 


5.70E+05 


9.80E+03 




R4 


6.84E-14 


8.29E+38 


6.06E+50 


-3.49 


1.66 


2.04E-03 


2.32E+05 


3.96E+03 




R5 


4.25E-14 


5.15E+38 


3.76E+50 


-3.65 


1.47 


1.37E-03 


4.71E+05 


2.12E+03 




R6 


6.43E-14 


7.79E+38 


5.70E+50 


-3.65 


1.81 


l.llE-03 


4.81E+05 


3.21E+03 




R7 


1.60E-13 


1.94E+39 


1.42E+51 


-3.58 


2.47 


8.31E-04 


7.11E+05 


7.03E+03 


NGC 3504 


R3+R4 


4.20E-13 


2.02E+40 


1.47E+52 


-3.32 


3.11 


6.76E-04 


4.70E+06 


8.07E+04 
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Figure 16. Relation between the N/O ratio and the O/H abundance for CNSFR (black circles and green upside down triangles for high 
and low metallicities respectively), high metaUicity HII regions (grey circles) and HII galaxies (cyan diamonds). 



with very little dependence on metaUicity or ionisation tem- 
perature. It can be seen that all the CNSFR observed show 
large [SII]/[SIII] ratios which imply extremely low ionisation 
parameters. On the other hand, a different answer would 
be found if the [OII]/[OIII] parameter, also commonly used 
as ionisation parameter indicator, was used. In this case, 
CNSFR and high metaUicity HII regions show a much more 
similar distribution, with CNSFR showing slightly lower val- 
ues of [OII]/[OIII] (Figure [Tt] right panel). 



Diaz et al. ( 1985) in order to investigate the possible contri- 



It should be noted the dependence of the [OII]/[OIII] 
parameter on metaUicity due to the presence of opacity edges 
of various abundant elements (O^, Ne^, C^^, N^^) in the 
stellar atmospheres that can combine to substantially mod- 
ify the stellar flux of high abundance stars at energies higher 
than 35-30 eV and then produce a lower [OIII] emission (Bal- 
ickfc Sneden|1976t . However, if the CNSFR were ionised by 
stars of a higher metaUicity than those in disc HII regions 
this effect would go in the direction of producing higher 
[OH] /[OIII] ratios for the CNSFR, and ionisation parame- 
ters derived from [OH] /[OIII] ratios would be found to be 
lower than those derived from [SII]/[SIII] ratios, contrary to 
what is actually observed. 

The ionisation structure can provide important infor- 
mation about the characteristics of the ionising source. A di- 
agram of the emission line ratios [OH] /[OIII] vs [SII]/[SIII], 
in particular, works as a diagnostics for the nature and tem- 
perature of the radiation field. This diagram was used in 



butions by shocks in CNSFR and LINERs and is the basis 
of the definition of the 77' parameter ( Vilchez Pagel|1988 ). 
The 77' parameter, defined as: 



[0//]AA3727, 29/[0///]AA4959, 5007 
[5'//] AA6716, 6731/[5'///]AA9069, 9532 



is a measure of the "softness" of the ionising radiation 
and increases with decreasing ionising temperature. The "77' 
plot" is shown in Figure [Ts] In this plot, diagonal lines of 
slope unity would show the locus of ionised regions with 
constant values of 77'. The lines shown in the plot have slope 
1.3 refiecting the second order dependence of 77' on ionisa- 
tion parameter ( Diaz et al.|1991| ). In this graph CNSFR are 
seen to segregate from disc HII regions. The former cluster 
around the value of logry' = 0.0 {Tion - 40,000 K) while the 
latter cluster around log 7y' = 0.7 {Tion ^ 35,000 K). Also 
shown are the data corresponding to HII galaxies. Indeed, 
CNSFR seem to share more the locus of the HII galaxies 
than that of disc HII regions. 

Given the problems involved in he measurement of the 
[OH] emission lines mentioned in section 4.2 it is interesting 
to see the effects that a possible underestimate of the inten- 
sity of this line would have on our analysis. The underesti- 
mate of the [OH] line leads to an overestimate of the SO 23 
parameter and hence to an underestimate of the [SHI] elec- 
tron temperature. This lower temperature in turn leads to 
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Figure 17. Distribution of the [SII]/[SIII] (left) and [OII]/[OIII] (right) ratios for the observed CNSFR (dark) and the sample of high 
metallicity disc HII regions (light). 
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Figure 18. The 77 ' plot: the [OII]/[OIII] ratio vs [SII]/[SIII] ratio for different ionised regions: high metallicity CNSFR (black circles), 
low metallicity CNSFR (green upside down triangles), high metallicity disc HII regions (grey circles and HII galaxies (cyan diamonds). 



higher 0/H abundances. Therefore, larger values of the in- 
tensity of the [Oil] line would lead to lower 0/H abundances. 
However, not by a large amount. Increasing the [OH] inten- 
sity by a factor of 2, would decrease the SO23 parameter by 
a factor between 1.5 and 1.8 depending on the region, which 
would lead to increased [SHI] electron temperatures by be- 
tween 500 and 1000 K. The corresponding 0^/H^ would be 
only slightly smaller, by about 0.10 dex, as would also be the 
total 0/H abundances which is dominated by 0^/H^. The 
0++/H+ ionic ratio however would be lower by between 0.3 
and 0.4 dex, which would produce an ionisation structure 
more similar to what is found in disc HII regions. The sul- 
phur ionic and total abundances would be decresed by about 
0.17 dex. Finally, the N+/0+ ratios would be decreased by 
about 0.17 dex, while the S/0 ratios would remain almost 
unchanged. In the "77' plot" (Fig. |18[) the data points corre- 
sponding to our observed CNSFR would move upwards by 
0.30 dex. However, we do not find any compelling reason 
why the [OH] intensities should be larger than measured by 
such a big amount (see section 4.2 and Fig. |8|. 

One possible concern about these CNSFR is that, given 
their proximity to the galactic nuclei, they could be af- 
fected by hard radiation coming from a low luminosity AGN. 
NGC 3551 shows a faint UV core . However, the lUE spec- 
trum that covers the whole central starforming ring, shows 



broad absorption lines of SilV A 1400 A and CIV A 1549 A 



typical of young stars of high metallicity ( Colina et al. 1997 ). 
They are consistent with a total mass of 3 x 10 M© of re- 
cently formed stars (4-5 Myr). This is of the order of our 
derived values for single CNSFR in this galaxy. Therefore, 
no signs of activity are found for this nucleus nor are they 
reported for the other two galaxy nuclei. On the other hand, 
the Hell A 4686 A line is measured in regions R1+R2 and 
R6 of NGC 2903 and in region R7 in NGC 3351. In the 
first region, there is some evidence for the presence of WR 
stars (Castellanos et al. 2002). For the other two, that pres- 
ence is difficult to assess due to the difficulty in placing the 
contiuum for which a detailed modelling of the stellar pop- 
ulation is needed. 

Alternatively, the spectra of these regions harbouring 
massive clusters of young stars might be affected by the 
presence of shocked gas. Diagnostic diagrams of the kind 



presented by Baldwin et al. (1981) can be used to investi- 
gate the possible contribution by either a hidden AGN or 
by the presence of shocks to the emission line spectra of 
the observed CNSFR. Figure ^] shows one of these diag- 
nostic diagrams, log ([Nil] /Ha) vs log ([0III])/H/3, for our 
CNSFR: light squares for NGC 3351, triangles for NGC 2903 
and dark square for NGC 3504. The figure has been adapted 



from Groves et al. (2006) and shows the location of emis- 
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sion line galaxies in the Sloan Digital Sky Survey (SDSS). 
Dashed and solid lines correspond to the boundary between 
Active Galactic Nuclei (AGN) and HII galaxies defined by 
Kewley et all (|2001| and Kauffmann et all ([2003) respec- 



tively. Some of our observed CNSFR are found close to the 
transition zone between HII region and LINER spectra but 
only R3+R4 in NGC 3504 may show a hint of a slight con- 
tamination by shocks. This region, the most luminous in our 
sample and also the one with the higest abundance, should 
be studied in more detail. 

A final remark concerns the gas kinematics in CNSFR. 
In a recent work l|Hagele et alT (2007) have studied the kine- 
matics of gas and stars in the CNSFR of NGC 3351 find- 
ing two different components for the ionised gas in H/3 and 
[OIII] emission: a "broad component" with a velocity disper- 
sion similar to that measured for the stars, and a "narrow 
component" with a dispersion lower than the stellar one by 
about 30kms~^. Obviously the abundance analysis and the 
location of these regions on diagnostic diagrams would be af- 
fected if more than one velocity component in the ionised gas 
corresponding to kinematically distinct systems are present. 



7 SUMMARY AND CONCLUSIONS 

We have obtained spectro-photometric observations in the 
optical (AA 3650-7000 A) and far-red (AA 8850-9650 A) 
wavelength ranges of 12 circumnuclear HII regions in the 
early type spiral galaxies: NGC 2903, NGC 3351 and 
NGC 3504. These regions were expected to be amongst 
the highest metallicity regions as corresponds to their po- 
sition near the galactic centre. At the same time, this posi- 
tion implies a substantial contribution by the bulge stellar 
population to their spectra which represents a major ob- 
servational problem and compromises the reliability of the 
emission line intensities. Its proper subtraction however re- 
quires the detailed modelling of the stellar population and 
the disentangling of the contribution by the bulge and by 
any previous stellar generations in the CNSFR themselves. 
Almost all these regions show the presence of the CaT lines 
in the far red (see Hagele et al. 2007) some of them with 
equivalent widths that suggest a certain contribution by red 
supergiant stars. A detailed study of the ionising and non 
ionising stellar populations of this regions will be presented 
in a forthcoming paper. In the case of the Balmer lines, the 
presence of wide wings in their profiles, allows to perform a 
two-component - emission and absorption - gaussian fit in 
order to correct the Balmer emission lines for this effect. 

We have derived the characteristics of the observed 
CNSFR in terms of size. Ha luminosities and ionising clus- 
ter masses. The derived sizes are between 1.5 arcsec and 
5.7 arcsec which correspond to linear dimensions between 
74 and 234 pc. The derived filling factors, between 6 x 10~^ 
and 1 X 10"*^, are lower than commonly found in giant HII 
regions (~ 0.01). Ha luminosities are larger than the typi- 
cal ones found for disc HII regions and overlap with those 
measured in HII galaxies. The region with the largest Ha lu- 
minosity is R3+R4 in NGC 3504, for which a value of 2.02 
X 10^° is measured. Ionising cluster masses range between 
1.1 X 10^ and 4.7 x 10^ M© but could be lower by factors 
between 1.5 and 15 if the contribution by the underlying 
stellar population is taken into account. 



The low excitation of the regions, as evidenced by the 
weakness of the [OIII] A 5007 A line, precludes the detection 
and measurement of the auroral [OIII] A 4363 A necessary 
for the derivation of the electron temperature. Only for one 
of the regions, the [SHI] A 6312 A line was detected pro- 
viding, together with the nebular [SIII] lines at A A 9069, 
9532 A , a value of the electron temperature of Te([SIII]) = 
8400^^250 A method for the derivation of sulphur 
abundances was developed based on the calibration of the 
[SIII] electron temperature vs the empirical parameter SO23 
defined as the quotient of the oxygen and sulphur abundance 
parameters O23 and S23 and the further assumption that 
T([SIII]) ~ T([SII]). Then the oxygen abundances and the 
N/0 and S/0 ratios can also be derived. 

The derived oxygen abundances are comparable to 
those found in high metallicity disc HII regions from direct 
measurements of electron temperatures and consistent with 
solar values within the errors. The region with the highest 
oxygen abundance is R3+R4 in NGC 3504, 12+log(0/H) = 
8.85, about 1.5 solar if the solar oxygen abundance is set at 
the value derived by Asplund et al. (2005), 12+log(O/H)0 
= 8.66±0.05. Region R7 in NGC 3351 has the lowest oxygen 
abundance of the sample, about 0.6 times solar. In all the 
observed CNSFR the 0/H abundance is dominated by the 
0^/H^ contribution, as is also the case for high metallicity 
disc HII regions. For our observed regions, however also the 
S^/S^^ ratio is larger than one, different from the case of 
high metallicity disc HII regions for which, in general, the 
sulphur abundances are dominated by S^^/H+. The derived 
N/0 ratios are in average larger than those found in high 
metallicity disc HII regions and they do not seem to follow 
the trend of N/0 vs 0/H which marks the secondary be- 
haviour of nitrogen. On the other hand, the S/0 ratios span 
a very narrow range between 0.6 and 0.8 of the solar value. 

When compared to high metallicity disc HII regions, 
CNSFR show values of the O23 and the N2 parameters 
whose distributions are shifted to lower and higher values 
respectively, hence, even though their derived oxygen and 
sulphur abundances are similar, higher values would in prin- 
ciple be obtained for the CNSFR if pure empirical methods 
were used to estimate abundances. CNSFR also show lower 
ionisation parameters than their disc counterparts, as de- 
rived from the [SII]/[SIII] ratio. Their ionisation structure 
also seems to be different with CNSFR showing radiation 
field properties more similar to HII galaxies than to disc high 
metallicity HII regions. The possible contamination of their 
spectra from hidden low luminosity AGN and/or shocks, 
as well as the probable presence of more than one velocity 
component in the ionised gas corresponding to kinematically 
distinct systems, should be further investigated. 
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Figure 19. The [NII]/Ha vs [OIII]/H/? diagnostics for emission line objects in the SDSS (adapted from [Groves et al.| p006t ). The 
location of our observed CNSFR is shown overplotted. Light (blue) triangles correspond to NGC 2903, light (yellow) squares correspond 
to NGC 3351 and the black square corresponds to NGG 3504. Solid ( [Kewley et ar]|2001| > and dashed ( [Kauffmann et al.|[2003t lines 
separate the starforming from the active nucleus galaxies regions. 
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